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2p-based magnetic moments and magnetic coupling are studied with density functional based 
methods for substitutional N in the alkaline earth monoxide series MgO, CaO, SrO, BaO. The 
hole is rather strongly localized near the N 2_ ion, being somewhat more so when strong on-site 
interactions are included in the calculations. Strong magnetic coupling is obtained in the itinerant 
electron limit but decreases strongly in the localized limit in which the Coulomb repulsion within 
the N 2p shell (U) is much greater than the N 2p impurity bandwidth (W). We find that this limit 
is appropriate for realistic N concentrations. Ordering on a simple cubic sublattice may maximize 
the magnetic coupling due to its high directionality. 

PACS numbers: 



I. INTRODUCTION 

After a strong research effort spent on magnetic- 
ion doped semiconductors as spintronics materials, 1 
scientific effort has recently turned to the possibil- 
ity of obtaining useful ferromagnetic (FM) order at 
room temperature in materials without conventional 
magnetic atoms with open d or / shells. Isolated de- 
fects in insulators often produce magnetic moments, 
which occur anytime a state in the band gap is oc- 
cupied by an odd number of electrons. A common 
example is the P impurity in Si, where the mag- 
netism of the material evolves with increased dop- 
ing as the insulator-metal transition is approached. 
Such strictly localized states in insulators are mag- 
netic regardless of their underlying atomic character 
or their degree of localization. 

Elfimov et ali£ proposed to use the cation vacancy 
- the charge conjugate of the conventional F-center 
state - in non-magnetic oxides such as CaO^ as a 
source of magnetic moment in insulators. The re- 
sulting pair of holes, which are bound to the region 
for an isolated vacancy, may form a magnetic cen- 
ter even if it does not lie in the gap, due to the 
large Hund's coupling on the oxygen ions neigh- 
boring the vacancy (which can no longer be simple 
non-magnetic O 2- ). These states extend to a few 
atomic neighbors, with magnetic coupling resulting 
from direct exchange. Magnetic order becomes pos- 
sible, and at concentrations of a few percent they 
suggested that conduction, and half metallic ferro- 
magnetism, might result. 

The Osaka group^ followed by suggesting that 
B, C, or N substituted for O in CaO could result 
in local moments in the 2p states of the impurity 
ions. The coupling between moments was found, in 
calculations on random alloys, to be FM and half- 
metallic ferromagnetism was also predicted, possibly 



as high as room temperature. Elfimov et al£ stud- 
ied the specific case of N substituting for O in SrO, 
applying not only the conventional local spin den- 
sity approximation (LSDA) but also the strongly- 
correlated form (LSDA+U) that should be appro- 
priate for localized magnetic states. They also ob- 
tained the possibility of half metallic ferromagnetism 
at technologically relevant temperatures. Experi- 
mentally, magnetism has been observed as a surface 
effect in thin films of these non-magnetic oxides^ 



TABLE I: Lattice parameters^ and band gap£ of the cor- 
responding monoxides under study. 



Compound 


Lattice parameter (A) 


Band gap (eV) 


MgO 


4.12 


7.2 


CaO 


4.81 


6.2 


SrO 


5.16 


5.3 


BaO 


5.52 


4.0 



There are other examples of atomic p charac- 
ter giving rise to magnetism in solids. One ex- 
ample is the class of alkali hyperoxides (for exam- 
ple, RbO^^ and Rb 4 0e±ii±2). These materials are 
ionic, containing the spin-half O^ - ion, and display 
a variety of unusual magnetism-related properties. 
The SrN2 compound, an ionic material containing 
the N 2- unit, is a related example?^ while the SrN 
compound with magnetic N is predicted to be a 
half metallic ferromagnet.— Between the bulk com- 
pounds and the (nearly) isolated magnetic impuri- 
ties, another promising example has been reported. 
At the p-type interface between LaA103 and SrTi03, 
there is (from simple electron counting) 0.5 too few 
electrons per interface cell to fill all of the O 2p 
states that are filled in both bulk materials. Holes 
mush exist in the O 2p bands at the interface, and 



Pcntcheva and Pickett predicted that magnetic holes 
arise and should order at low temperature. 15 Subse- 
quently, magnetic hysteresis has been observed at 
this interfaced 

Here we will analyze the magnetism and magnetic 
coupling in alkaline-earth monoxides doped with N 
in more detail. Substitution of O by N in an insulat- 
ing host presents a clear platform for local moments, 
and the basic picture is simple. Besides changing the 
potential locally, the substitution of O by N removes 
one electron giving a N~ 2 2p 5 ion, and thereby in- 
troduces a hole that is necessarily magnetic (being 
spin-uncompensated). Ferromagnetism may result 
if the magnetic coupling has the right character and 
is sufficiently robust. Within the series MgO, CaO, 
SrO and BaO the volume varies by a factor of 2.3 
and the bandgap varies by over a factor of 1.8, see 
Table HI By introducing dopant nitrogen atoms with 
concentrations of one N atom per 8, 16 and 32 O 
atoms, we are able to model impurity concentrations 
of 12.5%, 6% and 3%, respectively. The goal here is 
to understand the electronic structure and magnetic 
character of the substitutional N impurity, to study 
dopant interactions for the various concentrations, 
and to analyze the magnetic properties by calculat- 
ing the condition for strong coupling that could lead 
to a FM T c on the order of room temperature. 

These semiconductors crystallize in a rocksalt 
structure. Upon increasing the lattice parameter 
by introducing a bigger cation, there should be a 
tendency for the N hole state to become more local- 
ized. However, the band gap narrows in proceeding 
from the smallest to largest (7.2 eV for MgO, 4.0 for 
BaO), a trend that should cause a corresponding in- 
crease in the tendency towards delocalization of the 
impurity state. How these effects compete is impor- 
tant to determine. Another distinction is that the 
lowest conduction states in CaO, SrO, and BaO are 
d states, while MgO is very different in this regard 
(Mg s, O 3s, 3p). The hole state will be formed from 
valence orbitals, however, so the impact both of the 
magnitude of the band gap and the character of the 
conduction band must be calculated. 

Table |TT] summarizes the N-N nearest neighbor 
distances in the different compounds for the vari- 
ous concentrations considered. As the concentration 
changes, the coordination of the impurities varies, as 
is provided in the table. 



II. COMPUTATIONAL DETAILS 

Electronic structure calculations were performed 
within density functional theory-^ using wiEN2Kfi£ 
which utilizes an augmented plane wave plus local 
orbitals (APW+lo) 19 method to solve the Kohn- 



TABLE II: Nearest neighbors distances and number of 
nearest neighbors for the impurity at various concentra- 
tions and monoxides under study. 



Compound 


Concentration 


N-N distance (A) 
(no. neighbors) 


MgO 


1/8 


6.0 (12) 


MgO 


1/16 


6.0 (4) 


MgO 


1/32 


8.4 (6) 


CaO 


1/16 


6.8 (4) 


SrO 


1/8 


7.3 (12) 


SrO 


1/16 


7.3 (4) 


SrO 


1/32 


10.3 (6) 


BaO 


1/8 


7.8 (12) 


BaO 


1/16 


7.8 (4) 


BaO 


1/32 


11.0 (6) 



Sham equations. This method uses an all-electron, 
full-potential scheme that makes no shape ap- 
proximation to the potential or the electron den- 
sity. The exchange-correlation potential utilized was 
the Wu-Cohen version of the generalized gradient 
approximation^ and strong correlation effects were 
introduced by means of the LSDA+U scheme^ in- 
cluding an on-site effective U for the O and N p 
states. We have used option nldau= 1 in wien2k, 
i.e. the so-called "fully-localized limit" , using an ef- 
fective \] e ff— U-J, being J the on-site Hund's rule 
coupling constant, taken as J= 0. The justification 
of utilizing the LSDA+U scheme comes from spec- 
troscopic measurements in oxides, which estimate 
large values of U (5-7 eV) for the p states^ This 
has been shown to be crucial in understanding the 
interfacial electronic structure of transtion metal ox- 
ides, where magnetic moments can localize in the O 
atomsji 5 - All the calculations were converged with 
respect to all the parameters used, to the precision 
necessary to support our calculations. 



III. RESULTS 

We begin by studying the local electronic struc- 
ture of a N 2p 5 impurity in MgO. Comparing a non- 
magnetic solution for the N atom and a solution with 
1 /is/N (a local moment formation in the N site), 
the total energy difference favors the magnetic solu- 
tion by more than 100 meV/N, a value that is nearly 
independent of the concentration of impurities. The 
introduction of a value of U e // = 5 eV produces a 
further stabilization of the magnetic solution with 
respect to the non-magnetic one (the energy differ- 
ence rises above 300 meV/N). 
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FIG. 1: Total density of states for ferromagnetically 
aligned MgOi-^N^ with x= 1/16 ~ 6%, calculated with 
the LSDA+U method including spin-orbit effects using 
U= 5.5 eV for the 2p-states of both N and O. The ma- 
jority DOS is plotted upward, the minority is plotted 
downward. 
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FIG. 2: Total density of states of the various monox- 
ides under study for FM alignment, calculated with the 
LSDA+U method for U= 5.5 eV. Only the minority spin 
channel is shown, since it illustrates all the important 
features: O 2p bandwidth, narrowing of the bandgap, 
and the position of the hole state at 2-2.5 eV above the 
valence band maximum.. 



A. Effect of atomic relaxation 

Substitutional N in these monoxides will break the 
symmetry and may lead to some local distortion of 
the otherwise perfectly octahedral environment of 
the cation. To determine the importance of atomic 
relaxation we performed structural minimization for 
MgOi_ x N K , with x= 1/16. If sufficiently localized, 
the hole might occupy a p 2 orbital (1 2 = 0), leading to 
an elongation of the N-Mg distance along the z axis, 
or it could be p x ±ip y (1 2 = ± 1) that is unoccupied, 
leading to an elongation of the N-Mg distance in the 
plane perpendicular to the z axis. The symmetric 
LSDA solution shares the hole among the p orbitals 
equally. We have performed a structural optimiza- 
tion within GGA {i.e. without including strong cor- 
relation effects) which leads to an electron density 
describing the hole in the N atom being symmet- 
ric. The result is that Mg-N distances get modified 
only slightly, becoming elongated by less than 0.5%. 
Also, using the LSDA+U scheme for doing the struc- 
ture optimization, the hole located in a p 2 state leads 
to an elongation of about 3% of the Mg-N distance 
along the z-axis. We have found that these small 
differences have little effect on the trends we study, 
either the localization of the magnetic moment or 
the magnetic coupling, hence we have neglected re- 
laxation in the results that follow. 



B. Electronic and Magnetic Trends 

We have treated the N concentrations listed in Ta- 
ble|nl For each concentration, we have assumed that 



the N impurities are as regularly distributed as pos- 
sible (their separation is maximized) . As mentioned 
above, one can imagine how the physics can change 
from MgO to BaO as the size of the cation, and the 
volume, increases. On one hand, greater derealiza- 
tion is permitted as the band gap is reduced (by a 
factor 2 experimentally, by a factor 3 from our cal- 
culations for the doped compounds with U e //= 5.5 
eV). On the other hand, for a given concentration of 
impurities, the N-N distance becomes larger for the 
larger volume compounds. The former could encour- 
age an easier propagation of the magnetic coupling 
throughout the crystal, while the latter would im- 
ply a smaller magnetic interaction if the size of the 
defect state remains the same. 

We first consider the U e //= 5.5 eV electronic 
structure for a particular case, MgO with x= 1/16 ~ 
6%, where the density of states (DOS) for separate 
spin directions is shown in Fig. [TJ For the minority 
spin, the narrow hole band lies in the gap and has a 
full bandwidth of W « 0.3 eV. The other two N 2p 
minority states are occupied and lie above the O 2p 
bands. The separation of filled and unfilled minority 
N states is 2 eV, being the Hund's exchange about 
0.6 eV in the occupied N p states. The conduction 
bands show no appreciable spin splitting. When we 
include spin-orbit coupling, the hole occupies the p x 
- i p y state, l z = -1 orbital, satisfying Hund's rule. 

To display the trends across the series, it is suf- 
ficient to look only at the minority DOS, which is 
plotted for each compound in Fig. [2] The band 
gap reduction through the series Mg^Ba is evident. 
The O 2p bandwidth is much larger for MgO but 
this has little effect on the defect state, except per- 
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haps contributing a small additional broadening. In 
spite of this large band gap variation, and the change 
in the volume (hence the near neighbor distance) 
the position of the unoccupied hole state changes 
rather little, staying 2±0.5 eV above the highest oc- 
cupied state. The variation that does occur is non- 
monotonic and may be due to the conduction bands 
that fall in energy through the series. For BaO, nar- 
rowest band gap member, the N hole band merges 
into the bottom of the conduction Ba d bands. 

Because the intra-atomic Coulomb interaction 
strength U is at least a few eV, these systems are 
in the regime U/W >> 1 (being W the N 2p hole 
bandwidth) so hopping of the holes will be inhib- 
ited, giving insulating behavior. This is the limit in 
which the LSDA+U method works well. For values 
of U bigger than ~1 eV, the material is insulating 
whereas for smaller values (i.e. in LSDA or GGA), 
the N hole band intersects the upper valence bands 
and the band structure becomes metallic (or half- 
metallic). 

C. Sign, strength, and character of magnetic 
coupling 

Since the N ions form a strong local magnetic mo- 
ment, the next question is whether they couple fer- 
romagnetically and what the strength of the mag- 
netic interaction is. We have performed calcula- 
tions on both FM and antiferromagnetic (AF) align- 
ments to obtain the coupling between nearest neigh- 
bor N moments, for the different concentrations and 
monoxides under study. We estimate the Curie tem- 
perature by using the mean field expression from a 
Heisenberg model of the type H= - JijSiSj, with 
S=l/2 (the sum is over pairs), 

_ 2zS(S+l)J 

3k B ' [) 

where z is the number of nearest neighbors and J is 
the exchange constant. 

Table IIIII shows how the magnetic coupling be- 
tween N moments varies both with the concentra- 
tion and the different compound considered. The 
data presented in the table were obtained for a calcu- 
lation without considering strong correlation effects 
(U= 0), but including spin-orbit effects, and also 
using the Wu-Cohen GGA functional. The effect 
of introducing U is to reduce the magnetic coupling, 
because it further localizes the moments. In fact, us- 
ing a range of values of U up to 8 eV, we find that, for 
BaO and SrO, the introduction of U as small as 1.5 
eV leads to magnetic coupling too small to determine 
reliably from supercell energy differences. The hole 
states no longer overlap directly, and the materials 



TABLE III: Strength of the magnetic coupling converted 
into Curie temperature of the corresponding monoxides 
under study obtained from an LSDA calculation (includ- 
ing spin-orbit effects) and a GGA (Wu-Cohen) calcula- 
tion (U= 0). 



Compound 


Concentration 


Tc (K) 


Tc (K) 






[LSDA+SO] 


[GGA] 


MgO 


1/8 


330 


320 


MgO 


1/16 


400 


400 


MgO 


1/32 


5 


70 


CaO 


1/16 


72 


83 


SrO 


1/8 


82 


66 


SrO 


1/16 


31 


41 


SrO 


1/32 


250 


98 


BaO 


1/8 


32 


31 


BaO 


1/16 


11 


13 


BaO 


1/32 


280 


230 



are insulating with no means to transfer magnetic 
coupling except classical dipolar coupling. For MgO 
at the higher concentrations (x=l/8 and x= 1/16), 
the magnetic coupling is much higher than for the 
other compounds, because of the much smaller lat- 
tice parameter, and the introduction of U reduces 
the size of the interaction but does not kill T c . 

Putting all these LSDA data together, we can 
make a plot of how the magnetic interaction strength 
J varies in these compounds with respect to the N-N 
distance, neglecting all the differences between the 
various compounds we are considering. The result 
can be seen in Figure [31 that shows an exponential- 
like decay in the strength of the magnetic coupling 
as the N impurity atoms are separated. There is, re- 
markably, an upturn for N-N distances above 10 A, 
with values above room temperature for BaO and 
SrO for a concentration around x ~ 3%. For this 
supercell, the magnetic impurities are located on a 
simple cubic sublattice at lattice constant 2a. Elfi- 
mov et ai£ have shown that the electronic coupling 
between N impurities in SrO is highly directional, 
with strong coupling along the crystalline axes (large 
effective p x ~Px coupling along the x-axis, for exam- 
ple). This opens the possibility of reaching high tem- 
perature FM attainable concentrations of N dopant 
in the limit U~ W. 

A double-exchange-like microscopic mechanism 
seems to be the best explanation of the FM cou- 
pling in these monoxides in the GGA limit, when 
the Fermi level lies within the impurity band. Dou- 
ble exchange is normally associated with carriers 
coupled to a background of moments, with ki- 
netic energy favored by alignment of the moments. 
Electronic structure- wise, the added kinetic energy 
shows up as an increased bandwidth compared to 



4 



antialignment. Indeed for the GGA calculations and 
for small values of U, the system is (half) metallic 
and the double exchange picture as conventionally 
interpreted makes sense. As U is increased beyond 
~1 eV, a gap appears, and in the large U regime 
the coupling is certain to be via direct exchange. 
However, as shown in Fig. 3, the coupling strength 
(as reflected in T c in this figure) does not appear 
to undergo any discontinuity at the metal-insulator 
transition. The coupling (and T c ) does however de- 
crease rapidly as the value of U is increased. Having 
band states at the Fermi level leads also to a larger 
FM coupling. When U > U cr = 1 eV is introduced 
in the calculations, an insulating state is obtained 
and magnetic exchange coupling is reduced (see Fig. 

a. 



calculated the coupling for a variety of concentra- 
tions, paying particular attention to the effect of 




U (eV) 




FIG. 3: Exchange coupling strength with respect to the 
N-N distance, considering all the LSDA data from the 
various compounds together. "Error" bars simply denote 
the range of values that were calculated for the same near 
neighbor distance, see Table II. 



In Fig. 2] we illustrate more explicitly how the 
effect of U reduces the Curie temperature. Strong 
correlation effects produce a stronger localization of 
the magnetic moments, i.e. their effects extend to 
fewer neighbors. The reduced separation of overlap 
leads to a corresponding reduction of the magnetic 
interaction. Hence, in the limit U >> W, where this 
system resides for the concentrations we study, the 
magnetic coupling will not be strong. 



IV. SUMMARY 

In this paper we have studied the formation of lo- 
cal moments when substitutional N impurities are 
introduced in alkaline-earth monoxides. We have 



FIG. 4: Strength of the magnetic coupling in MgO ob- 
tained from an LSDA+U calculation (including spin- 
orbit coupling) with various values of U calculated in 
terms of the Curie temperature. 



Hubbard repulsion U (since U > W for these sys- 
tems, where W is the impurity band width). The 
local moments that form, mainly confined to the N 
atoms and their nearest neighbors, would lead to 
room temperature ferromagnetism in the itinerant 
electron limit via a double-exchange-like exchange 
coupling within the impurity band. However, in 
the more realistic localized electron limit U >> W 
the coupling must be due to direct exchange (lo- 
cal impurity state overlap) and becomes drastically 
reduced. Interestingly, the crossover between the 
double-exchange coupling and the direct exchange 
as the interaction strength is increased does not 
lead to any discernible anomaly in magnetic cou- 
pling strength at the metal-insulator transition. Fi- 
nally we confirm that the impurity state is highly 
anisotropic, and as result magnetic coupling is en- 
hanced by ordering of the magnetic N impurities on 
a periodic (simple cubic) sublattice. 
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